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THE CONNECTION BETWEEN CONVECTIVE DOWNFLOWS AND SPINES IN A SUNSPOT PENUMBRA 
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ABSTRACT 

We discuss NICOLE inversions of Fe I 630.15 nm and 630.25 nm Stokes spectra from a sunspot penumbra, 
recorded with the CRISP imaging spectropolarimeter at the Swedish 1-m Solar Telescope and a spatial resolu- 
tion close to 0715. Our emphasis is on narrow downflow lanes, which are cospatial with the relatively dark and 
cool parts of penumbral filaments. We find that these downflow lanes are located at the boundaries between 
areas of relatively horizontal magnetic field (the inter-spines) and much more vertical field (the spines). These 
locations agree with predictions from the convective gap model (the "gappy penumbra") proposed six years 
ago, and more recent 3D MHD simulations. We also find that in the deep photosphere, the downflow lanes are 
frequently, but not always, associated with opposite polarity magnetic field. This is also consistent with the 
simulations, predicting that many of the convective downflows pull down the magnetic field, rather than flow 
passively along it. 

Subject headings: Sunspots — Convection — Sun: magnetic topology — Sun: photosphere — Sun: surface 
magnetism 



1. INTRODUCTION 

The strong, nearly horizontal and radial outflo ws discov- 
ered more than 100 years ago dEvershedl 119091) . have in- 
spired 1-dimensional (ID) numerical modeling of the penum- 
bral filaments as flux tubes, with a siphon flow indi- 
rectly driven by the assumed difference in magnetic field 
strength between the ascending (in the penumbra) and de- 
scendi ng (outside the penumbra) parts of magnetic flux 
tubes dMever & Schmidull968HMontesinos & Thomasll 19891 
1 1 9971) . ISchlichenmaier et all d 1 998b! la|) instead modelled these 
flows as transient flows triggered by heating of thin flux tubes 
embedded in a convectively unstable penumbral atmosphere 
of given properties. However, neither of these models can 
provide the overall heat fl ux needed to explain the radia- 
tive losses of the penumb ra (ISchlichenmaier & Solankill2003l : 
ISpruit & Scharme^ 120061) . For this, efficient convection is 
needed, with vertical RMS velocities comparable to those of 
th e quiet Sun. 

ISpru it & Scha rmerl d2006l) proposed that the penumbral 
filamentary structure and complex magnetic field topology 
is the result of convection opening up radia l ly ali gned, 
essentially field-free gaps. ISanchez Almeidal d2005l) and 
ISanchez Almeida & Ichimotol ( 120091) proposed to explain ob- 
served penumbral Stokes profiles in terms of spatially unre- 
solved fluctuations in velocity and magnetic field, the Micro- 
Structu red Magnetic Atmosphe r es (MlSMAs). Recent simu- 
lations dHeinemann et al.l 120071: IRempel et ail 120091: iRempell 
1201 ll 120 121) clearly support the presence of strong penum- 
bral convection at observable scales, and lead to the sur- 
prising conclusion that the Evershed flow is identical to 
the ho r izontal component o f this convection dScharmer et alJ 
l2008ah IRempel et al.ll2009t) . This interpretation of penum- 
bral fine structure is however not without controversy, the 
main argument being that observational evidence for con- 
vective downflows we ll inside the outer boundary of the 
penumbra is missing dFranz & Schlichenmaierl 12009). Re- 
cently, such evidence was reported in the C I 538.03 nm 



line dScharmer etaTI 1201 U Uoshi et alJ 1201 ll) an d the Fe I 
630.15 nm line dScharmer & Henriquesl 120121) . based on 
observations w i th the Swedish 1-m Solar Telescope (SST; 
IScharmer et alJ 120031) and its imaging spectro polarimeter 
CRISP dScharmer||2006t IScharmer et"atll2008bl). Synthetic 
spectra calculated from simulations dRempe]||201 ll) confirm 
that such downflows should be visible in the C I 538.03 nm 
line at the spatial resolution of the SST dBharti et all 120 111) . 
In the present paper, we apply Stokes inversions to Fe I 
630.15 nm and 630.25 nm spectra to establish the morpho- 
logical connection between the convective downflows and the 
magnetic field in a sunspot penumbra. 

2. OBSERVATIONS AND DATA REDUCTION 

The observed data analyzed in this Letter are of a rea- 
sonably symmetric sunspot at approximately 15° heliocen- 
tric distance, obtained on 23 May 2010 with the SST 
and its imaging spectropolarimeter C RISP. These data are 
identical to those described recently dScharmer et alJ |201 It 
IScharmer & Henriqu es 2012), but here we include the 630.2 
nm line in the analysis and a continuum wavelength 60 pm 
redward of the 630.2 nm line, adding to the 15 wavelengths 
scanned for each of the Fe I 630.15 and 630.25 nm lines. A 
complete scan of the two Fe I lines required only 16.5 s ob- 
serving time. 

We first compensated the imag es for darks a n d the "raw" 
flats similar to as described by ISchnerr et alJ (1201 ll) . and 
then processed the im ages with the MOMFBD code of 
van Noort et all ( 120051) . This was followed by a small-scale 
dewarping of the images to remove remaining alignment er- 
rors from high-altitude seeing at arc second scale, which can- 
not be compensated for by the 4" x 4" subfield MOMFBD 
processing (Henriques 2012). We derived the transmission 
profile of CRISP and the pre-filter at each pixel, as well as 
a "true" flat-field image at each wavelength, by comparing 
the observed flat-field spectra, obtained close to disk cen- 
ter, to the FT S atlas spectrum (Brault & Neckel 1987; see 
lNeckellll999l) . The method used goes beyond that described 
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Figure 1. The left column (top to bottom) shows the temperature (scaled from 5000 K to 6600 K) at r c = 1, the LOS velocity obtained at r c = 1 from inversions 
applied to the Fe I 630. 15 nm and 630.25 nm lines, and the line core velocity of the C 1 538.03 nm line. The right column shows (top to bottom) the LOS magnetic 
field (saturated at -1000 G and 1400 G) obtained at r f = 0.02 and r c = 1 from inversions of the Fe I lines, and the signed integrated circular polarization of the 
C I 538.03 nm line. Tick marks are at 1" intervals, white dots are separated by 2". The arrow points in the direction of disk center. The FOV is ~ 35 X 19 ". 



in lScharmer et al.l d2008bl) by using analytical expr essions for 
the tr ansmission profiles of the two etalons (c.f., IS charmer! 
120061) . including the first side lobe on each side of the main 
peak for the high-resolution (HR) etalon. The free parameters 
of the fit are the wavelength shifts of the two etalons, the re- 
flectivity of the HR etalon and the parameters of an assumed 
quadratic variation of the pre-filter transmission with wave- 
length. Including the side lobes of the HR etalon transmis- 
sion profile returned excellent fits and an average reflectivity 
(93.44%) that is very close to that given by the manufacturer 
of the etalon coating (93.56%), suggesting very low levels of 
spectral straylight. From the fits, we also determine and com- 
pensate for the telluric blend in the red wing of the 630.25 nm 
line. In particular, this allows the inversions to take into ac- 
count any asymmetries in the transmission profile resulting 
from the relative wavelength shifts of the peak transmissions 
of the two Fabry-Perot etalons. Such asymmetries could oth- 
erwise lead to errors in the LOS velocity gradient. 

We also re-determined the polarization properties of the 
SST, using calibration data obtained with a 1-m rotating lin- 
ear polarizer mounted in front of the SST in May, June and 



October 201 1 . After demodulation of the MOMFBD restored 
images with an improved version of the polarization model of 
LSelbind (12001 . we checked the data for residual cross-talks 
from / to Q, U and V and from V to Q and U using meth- 
ods sim ilar to those described bv lSchlichenmaier & Colladosl 
( 120021) . The cross-talks found were 0.3% or less from / to Q, 
U and V, 2.7% from V to U, and much smaller from V to Q. 
These cross-talks were compensated for. 

As the final step of our pre- processing, we com- 
pensated the data for stra y-light (IScharmer et al.l 1201 It 
IScharmer & Henriquesll20f2l) . We assume the following re- 
lation between the observed I„ and "true" /, intensities at any 
wavelength and polarization state: I = (1 - a)I t + al t * G(W), 
where a is the straylight fraction, " * " denotes convolution, 
and G is a Gaussian straylight point-spread function (PSF), 
having a full width at half maximum of W. For the present in- 
versions, we set W to 1'.'8 and a to 0. 4, rather than to l'/2 and 
0.5 re spectively, as used previously ( IScharmer & HenriquesI 
2012]). Our present choice of straylight parameters give sim- 
ilar RMS contrast for granulation, and similar up- and down- 
flow velocities in the penumbra, but result in fewer failures 
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with the inversion code. The wavelength scale (used for es- 
tablishing the LOS velocity reference) was set by compar- 
ing the observed line profiles for granulation surrounding the 
sunsp ot, and void of strong fields, with those of 3D simula- 
tions (Ide la Cruz Rodrigue z et al. 2 01 ll) . 

3. INVERSIONS 

Inversions were applied jointly to the two Fe I lines. Abun- 
dances and atomic pa rameters used are identical to those of 
ISocas-Navarrol(l201 ll) . Inversions were made assuming a sin- 
gle atmospheric component in each pixel, excluding straylight 
compensation beyond that of the pre-processing (Sect. 2). 
We used the inversion code NICOLE, which is an improved 
implem entation of the inversion code of ISocas-Navarro et alj 
( 120001) . NICOLE was modified to accept a separate transmis- 
sion profile for each pixel and was applied iteratively to the 
data with an increasing number of free param eters (FP), sim- 
ilar to as described by ISocas-Navarrol (1201 ll) . We used the 
following final number of FP: 3 each for temperature and the 
LOS velocity, 2 each for B x , B v and B z , and 1 for the micro- 
turbulence (13 FP). We used linear/quadratic interpolation in 
log t c of the corresponding atmospheric variables with 2/3 FP. 
We also made inversions with 3 FP for B z in order to check 
that the inferred polarity reversals (Sect. 5) at continuum op- 
tical depth t c =1, are robust results, which is the case. Indeed, 
we find that the opposite polarity patches seen with 2 FP for 
B z are stronger and/or spatially more extended when using 3 
FP. We also see a few more opposite polarity patches in the 
inner penumbra with 3 than with 2 FP. However, the inver- 
sions with 3 FP for B z show enhanced levels of random noise 
and we discuss here only the results obtained with 2 FP for 
B z . We note that the fina l number of FP used in our inversions 
is smaller than used bv ISocas-Navarrol d201 ll) . This is a di- 
rect consequence of the spectral resolution of CRISP (about 
6 pm at 630 nm) being smaller than for his SOT/Hinode data, 
reducing the "height" resolution along the LOS for our data. 
We made a total of 41 inversion experiments in order to ex- 
plore the effects of changes in the pre-processing of the data 
and the number of FP, and spent considerable time verifying 
the robustness of our results. 

Figure[T]shows the LOS component of the magnetic field in 
the observers frame. Using the angle between solar North and 
disk center and the heliocentric distance of the spot, we trans- 
formed the magnetic field to the local frame. We resolved the 
180 degree ambiguity by defining an approximate center for 
the sunspot and chosing the azimuth angle offset (0 or 180 
deg) for which the horizontal field at r c =0.1 is directed more 
outward than inward from the spot center (this procedure for 
resolving the 180 degree ambiguity is of course inappropriate 
outside the sunspot). Figure [2] shows the magnitude and ori- 
entation of the horizontal field for every tenth pixel, and the 
vertical field (color coded) in the local frame. As is evident, 
the horizontal field is approximately parallell to the penum- 
bral filaments and there are no indications of major systematic 
errors from the above described procedure, or from the cali- 
bration of and compensation for the telescope polarization. 

4. GENERAL APPEARANCE OF THE INVERSION RESULTS 

As far as we know. lBellot Rubio et alj (120061) made the first 
measurements of LOS velocity (and temperature) gradients 
in a sunspot penumbra using inversion techniques applied to 
data from an imaging spectrometer. Here, we estimate gra- 
dients of both the magnetic field and the LOS velocity from 
CRISP data for the first time, and the question arises whether 



the results are meaningful and credible. We have investigated 
the velocity and magnetic maps returned by NICOLE over the 
entire 55" x 55" FOV and note the following: 

1) Isolated magnetic bright points and the interiors of larger 
flux concentrations (plage and pores) show a LOS magnetic 
field that decreases in strength with height, whereas at the 
boundaries of plage field the LOS magnetic field increases in 
strength with height. These results are consistent with obser- 
vations, models and simulations of flux tubes expanding with 
height (e.g..lSolanki et al.lll999tlYelles Chaouche et alJl2009b 
iPietarila et al.ll20 lOL and references therein). 

2) In the umbra, the light bridge and several bright umbral 
dots show a LOS magnetic fi eld increasing in strengt h with 
height (as earlier described by iRiethmuHer et al.ll2008l) . con- 
sistent with a field that wraps around and closes above struc- 
tures that have weak field, or are nearly field-free, deeper 
down (c.f., Figs.[Tp andQJl for examples). 

3) In the penu mbra, the interiors of magnetic spines 
dLites et alJI 19931) show a LOS magnetic field that decreases 
in strength with height (as is evident by comparing Figs. QJ> 
andQJ), but at the boundaries of the spines, the LOS magnetic 
field strength increases with height, consistent with canopy 
fields folding over the inter-spines (the "gaps"). Similar mag- 
netic f ield topologies were first described by iBorrero et al.l 
d2008h . 

4) Examination of the map of the transverse field (Fig. |2]i 
shows a clear tendency for the transverse field to diverge away 
from the penumbral spines, qualitatively consistent with flux 
conservation and the strong decrease of the strength of the 
LOS field with height in the spines. 

5) For the brightest umbral dots, the transverse field is much 
stronger on the disk center side than on the limb side in the 
deep layers, consistent with a field that folds over these dots. 

6) The measured LOS velocities decrease in strength with 
height al most everywhere in the FOV, in qualitative agreement 
with e.g lFranz & Schlichenmaieri (| 2009() . and as expected for 
overshooting convection. 

We will discuss some of these aspects of our Fe I data in forth- 
coming publications. 

Figure [T] shows the 35" x 19" sub-field discussed in this 
Letter. This Figure illustrates that the Fe I LOS velocity field 
at t c — 1 (Fig. Ub) in the quiet Sun and penumbra is very 
similar to that of the C I line core (Fig. U)p), form ed in the 
deep photosphere dSchlichenmaier & Schmidt! 19991) . We also 
note that the penumbral LOS magnetic field maps at t c = 1 
(Fig.\l)d) and the deeply formed C I line (Fig.\Tf) show a num- 
ber of nearly co-spatial opposite polarity patches, whereas in 
many cases no such opposite patches can be seen in the Fe I 
lines at t c = 0.02 (Fig.QJi). We thus find several independent 
pieces of evidence supporting our claim that the inversions re- 
turn trustworthy information about the height variations of the 
LOS velocity and magnetic field. 

5. DARK DOWNFLOW LANES 

Figure [3^ shows the temperature map at t c = 1, with the 
narrow dark redshifted lanes seen at r c = 1 in Fig.QJ overlaid 
in red. Several such lanes are clearly visible where the radial 
direction of the filaments is close to +90° from the disk center 
direction (the uppermost parts of the panels). At these loca- 
tions, there are no contributions from the radial Evershed flow, 
such that we conclude that these narrow and long redshifted 
lanes represent downflow lanes. Many similar looking narrow 
lanes, but with stronger redshifts, can be seen at the limb side 
penumbra, strongly suggesting that they are also downflow 
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Figure 2. The Figure shows the orientation (direction of the arrows) and strength (proportional to the lengths of the arrows) of the horizontal magnetic field at an 
optical depth of t c = 0.1, after transformation to the local frame. The vertical magnetic field (also in the local frame) is shown at r c = 1 and is encoded in color 
such that red corresponds to field of the same polarity as the spot and stronger than about 1500 G, and violet to opposite polarity field stronger than about 300 G. 
The arrows are plotted for every tenth pixel (corresponding to 0759), and give only a low-resolution view of the horizontal field. Tick marks are at 1" intervals. 
The FOV is ~ 35 X 19 " and corresponds to that shown in Fig.[T] 




Figure 3. The top panel shows the temperature (scaled from 5000 K to 6600 K) at r c = 1, with the dark downflow lanes at t c = 1 in Fig. 1, outlined in red. 
The second panel shows the same downflow lanes overlaid on a map of the vertical component (in the local frame) of the magnetic field (scaled as in Fig. 1) at 
t c = 1 . The third panel shows the vertical magnetic field, with opposite polarity field less than -200 G indicated in green. The fourth panel shows the magnitude 
of the horizontal field (in the local frame) at r c = 1, scaled from 150 G to 1800 G. 
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Figure 4. Observed (dashed) and fitted synthetic (solid) Stokes profiles in a 
downflow lane pixel, and the corresponding variations of temperature, LOS 
velocity, vertical (solid) and horizontal (dashed) magnetic field in the local 
frame). The 630.15/630.25 nm line profiles (plotted side by side) are normal- 
ized to the average continuum intensity outside the spot, and are shown for 
the pixel marked with "+" in Figs.fJJand [3] 

lanes and that they are associated with strong outflows. This 
conclusion receives strong support by the absence of narrow 
blueshifted flow lanes on the disk center side. If the narrow 
lanes on the limb side instead were strong outflow lanes with 
(weak) upflows, they would show stronger Doppler signature 
(but with opposite sign) on the disk center side, where the 
LOS would be more or less aligned with the flow direction, 
than on the limb side. Since the opposite is true, we feel con- 
fident in our identification of these lanes as downflow lanes. 

Figure|4]shows the observed and fitted synthetic Stokes pro- 
files for a downflow lane pixel marked with "+" in Figs. Q] 
and [3] The Stokes V profiles are abnormal with 3 lobes in- 
stead of 2, but are nevertheless well fitted by the variations 
of the LOS velocity and vertical magnetic field (in the local 
frame) with t c , shown in the lower panels. The polarity of the 
vertical magnetic field reverses below r e = 0.2, and the mag- 
nitude of the horizontal field decreases rather strongly with 
height. 

Comparing the temperature maps in Figs. QJ and [3^, it is 
clear that almost all downflow lanes well inside the outer 
boundary of the penumbra are co-spatial with the relatively 
(compared to their immediate surroundings) cool parts of 
the penumbral filaments, blinking the two maps makes this 
even more obvious (the opposite is not true: there are many 
relatively cool filamentary structures that are not associated 
with obvious downflow lanes). We conclude that these down- 
flow lanes are part of the convective flows, and contribute to 
the intensity-v ertical velocity correlation established in our 
previo us work dScharmer et al. 201 it IScharmer & Henriquesl 
l2012h . 

Figure[3]3 shows a remarkable relation between the location 
of the downflow lanes and the magnetic spines (shown white): 
well inside the penumbra, most of the downflow lanes are lo- 
cated at the boundaries of the spines, or, where the spines are 
very close together, between the spines. 

Figure [3J; shows in green color locations where the verti- 
cal magnetic field component (in the local frame) is of op- 
posite polarity with a strength of more than 200 G. Most of 
this opposite polarity field is located close to the outer bound- 
ary of the penumbra and is associated with (in some cases 
very strong) downflows. In the interior parts of the penumbra, 
many patches of opposite polarity field are co-spatial with the 



narrow downflow lanes. 

The magnitude of the horizontal field at t c =1 shown in 
Fig-EJl demonstrates that the spine structure, which is obvious 
in the map of the vertical field, is absent in the horizontal field 
map. This means that the horizontal field is nearly as strong 
in the spines as in the inter-spines at this depth. 

6. DISCUSSION 

We have identified narrow downflow lanes in the interior 
penumbra that are co-spatial with the relatively cool parts of 
the penumbral filaments, and located at the boundaries of the 
magnetic spines. These locations agree with predi ctions of the 
convective g ap model (IScharmer & Spruitj 120061) and recent 
simulations dRempell l 2012l his Fig. 22). We also find that in 
most cases these downflow lanes are associated with opposite 
polarity field, supporting the identification of these lanes as 
downflow lanes. The absence of opposite polarity field asso- 
ciated with some downflow lanes also agrees with the simu- 
lations, showing that the downflows frequently pull down the 
magnetic field instead of just flowing passively along it. 

The horizontal field inferred from our data is nearly as 
strong in the inter-spine s as in the spin es. This is consistent 
with the simulations of iRempel d2012l) . showing a layer of 
strong horizontal magnetic field near r c = 1, and weak-field 
gaps located below the visible surface. The distinctly different 
morphologies of the penumbral horizontal and vertical field 
components show n in Fig. [3] are remarkably similar to those 
of the simulations (Rempel 2012, his Figs. 6b and 6c). 

By calibrating the CRISP transmission profile at each pixel, 
we have been successful in using the inversion code NICOLE 
to simultaneously estimate LOS velocities, the magnetic field 
vector, and their LOS gradients from combined 630.15 nm 
and 630.25 nm Stokes spectra, at a spatial resolution close 
to 0'.'15. These results support earlier measurements in the 
deeply formed C I line. The present and earlier discussed 
(IScharmer et alj|20lll : IScharmer & Henriquesl2012l) observa- 
tional data, as well as simulations (lRempelll2012l). and svn - 
thetic spectra obtained from simulations dBharti et alj|201 ll) . 
clearly illustrate the need for reaching this spatial resolution 
in order to resolve the convective downflows in the interior of 
sunspot penumbrae. 
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